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Cj ' Abstract. 

r ■) , We present a new model for the infrared emission of the high redshift hyperluminous infrared 

. ' galaxy IRAS F10214+4724 which takes into account recent photometric data from Spitzer and 

iS^ ' Herschel that sample the peak of its spectral energy distribution. We first demonstrate that the 

t^H, combination of the AGN tapered disc and starburst models of Efstathiou and coworkers, while 

Q able to give an excellent fit to the average spectrum of type 2 AGN measured by Spitzer, fails to 

$H , match the spectral energy distribution of IRAS F10214+4724. This is mainly due to the fact that 

^ • the vSv distribution of the galaxy falls very steeply with increasing frequency (a characteristic 

C^ ' of heavy absorption by dust) but shows a silicate feature in emission. We propose a model that 

assumes two components of emission: clouds that are associated with the narrow-line region 

. • and a highly obscured starburst. The emission from the clouds must suffer significantly stronger 

►^ ' gravitational lensing compared to the emission from the torus to explain the observed spectral 

l-*^ , energy distribution. 

SJ^ . Keywords, galaxies: evolution, galaxies: formation, galaxies: high-redshift, infrared: galaxies, 
^-.^ ' cosmology: observations 

rn ! 

CN ' 1. Introduction 

T! ■ The title of the discovery paper of IRAS F10214+4724 at a redshift of 2.285 (Rowan- 

K^ ' Robinson et ah 1991) was 'An IRAS galaxy virith huge luminosity - hidden quasar or 
protogalaxy?'. It soon became clear that the reason for the enormous luminosity was 
5^ that the galaxy was gravitationally lensed by a factor of 10-100 (e.g. Broadhurst & Lehar 

Ci • 1995). Until recently, however, it was not clear whether the dominant energy source was 
an active galactic nucleus (AGN) or a starburst. 

Models for the infrared spectral energy distribution (SED) of IRAS F10214-f4724 in 
the mid-90s remained inconclusive (e.g. Rowan-Robinson et al. 1993) mainly due to 
the limited wavelength coverage. Mid-infrared spectrophotometry has been shown for 
some time now to be able to differentiate between a starburst and AGN origin of the 
infrared emission (Roche et al. 1991 and references therein, Genzel et al. 1998). The 
mid-infrared spectrum of IRAS F10214+4724 obtained by Teplitz et al. (2006) using the 
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IRS onboard Spitzer showed a puzzUng emission feature at rest frame lO/irn instead of 
the expected absorption feature given the type 2 nature of the AGN suggested by optical 
observations. Efstathiou (2006) proposed that the SED of IRAS F10214+4724 could be 
understood in terms of three components of emission: clouds that are associated with 
the narrow- line region and with a covering factor of 17%, an edge-on torus and a highly 
obscured starburst. It was also noted by Efstathiou (2006), however, that because the 
emission from the torus and the narrow-line region clouds arises from diferent regions 
their magnification may be different. 

In this paper we develop a new model for the infrared emission of IRAS F10214-I-4724 
which takes into account new photometric data from Spitzer and Herschel that cover the 
wavelength range 70-170 iJ.m and therefore complete its SED. 



2. Radiative transfer models 

2.1. Starburst models 

Efstathiou, Rowan-Robinson & Siebenmorgen (2000) developed a starburst model that 
has three main features: The first feature of the model is that it incorporates the stellar 
population synthesis model of Bruzual & Chariot that gives the spectrum of the stars 
as a function of their age. Efstathiou, Rowan- Robinson & Siebenmorgen (2000) used the 
table that assumes a Salpeter IMF and stellar masses in the range 0.1-125 Mq. The 
second important feature of the model is that it involves detailed radiative transfer that 
takes into account multiple scattering from grains and incorporates a dust model that 
includes small transiently heated grains and PAH molecules as well as large classical 
grains. For the calculation of the emission of the small grains and PAHs the method of 
Siebenmorgen & Kriigel (1992) was employed. The third feature of the model is that it 
incorporates a simple model for the evolution of the giant molecular clouds that constitute 
the starburst once a stellar cluster forms instantaneously at the center of the cloud. The 
most important characteristic of this evolutionary scheme is that by about 10 Myrs 
after star formation, the expansion of the HII region leads to the formation of a cold 
narrow shell of gas and dust. This naturally explains why the mid- infrared spectra of 
starburst galaxies are dominated by the PAH emission and not by the emission of hot 
dust. The sequence of molecular clouds at different ages can be convolved with a star 
formation history to give the spectrum of a starburst at different ages. Farrah et al. 
(2003) showed that the starburst models of Efstathiou, Rowan-Robinson & Siebenmorgen 
(2000) in combination with the tapered discs of Efstathiou & Rowan-Robinson (1995), 
which are described below, provided good fits to the SEDs of 41 ultraluminous infrared 
galaxies. One interesting prediction of the Efstathiou, Rowan-Robinson & Siebenmorgen 
(2000) model is that young starbursts show strong near-infrared continuum, small 6.2/i?Ti 
equivalent widths and deep silicate absorption features and can therefore explain the 
galaxies in class 3A of Spoon et al. (2007). The starburst model has recently been revised 
by Efstathiou & Siebenmorgen (2009) 

2.2. AGN torus models 

The most important constraint on early models for the torus in AGN (Pier & Krolik 
1992, Granato & Danese 1994, Efstathiou & Rowan-Robinson 1995) was provided by mid- 
infrared spectroscopic observations from the ground in the 8-13fim window by Roche, 
Aitken and collaborators (Roche et al. 1991 and references therein). The observations 
showed moderate absorption features in type 2 AGN and featureless spectra in type 1 
AGN. By comparison the flared discs (whose thickness increases linearly with distance 
from the central source) of Efstathiou & Rowan-Robinson (1995)) and Granato & Danese 
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(1994) showed strong silicate emission features when observed face-on. Efstathiou & 
Rowan-Robinson (1995)) proposed that tapered discs (whose thickness increases with 
distance from the central source in the inner part of the disc but tapers off to a constant 
value in the outer disc) with a density distribution that followed r~^ could give flat spectra 
in the mid-infrared for an opening angle of around 45° and an equatorial optical depth 
at lOOOA of 1000. Models with a smaller opening angle predicted face-on spectra with 
shallow silicate absorption features whereas models with a larger opening angle showed 
emission features. Recent Spitzer observations showed emission features in quasars (e.g. 
Hao et al. 2005, Hao et al. 2007, Siebenmorgen et al 2005, Spoon et al. 2007). The tapered 
discs of Efstathiou & Rowan-Robinson (1995), in combination with the starburst models 
of Efstathiou, Rowan- Robinson & Siebenmorgen (2000) have been very successful in 
fitting the spectral energy distributions of a number of AGN even in cases where mid- 
infrared spectroscopy was available (Alexander et al. 1999, Ruiz et al. 2001, Farrah et al. 
2003, Efstathiou & Siebenmorgen 2005). 

In this paper we use a grid of tapered disc models computed with the method of 
Efstathiou & Rowan- Robinson (1995). In this grid of models we consider four discrete 
values for the equatorial lOOOA optical depth (500, 750, 1000, 1250), three values for 
the ratio of outer to inner disc radii (20, 60, 100) and three values for the opening angle 
of the disc (30°, 45° and 60°). The spectra are computed for 37 or 74 orientations which 
are equally spaced in the range to tt/2. 

The grid of tapered disc models, as well as the grid of starburst models discussed 
above, are available from AE on request (a.efstathiou@euc.ac.cy). 

2.3. Comparison with the average spectrum of Seyfert 2s 

We first explore whether the combination of the tapered disc and starburst models de- 
scribed above can fit the average spectrum of Seyfert 2s presented by Hao et al. 2007. 
The observed spectrum has been determined by first normalizing the spectra at 15/im 
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Figure 1. Fit to the spectral energy distribution of IRAS F10214-I-4724. The dotted and dashed 
lines show the contributions of the AGN and starburst respectively to the total emission (solid 
line). The dash dotted line also shows a fit to the IRAC data with a IGyr old stellar population. 
Data from Teplitz et al. (2006), IRAS, Rowan-Robinson et al. (1993), Sturm et al. (2010). 



4 Andreas Efstathiou, Natalie Christopher, Aprajita Verma & Ralf Siebcnmorgen 

and then averaging. The AGN models have also been averaged in an attempt to mimic 
the averaging of the observed spectra. 

We find that the best fit to the Seyfert 2 average is provided by a model that assumes 
a torus opening angle of 45°, an equatorial lOOOA optical depth of 750 and a ratio of 
outer to inner torus radii of 20. The best fit starburst model assumes an initial optical 
depth of the molecular clouds Ty of 200 and a starburst age of 72Myr. 



3. A new model for IRAS F10214+4724 

We find that the models described above provide a very poor fit to the SED of IRAS 
F10214+4724. This is clearly due to the fact that the SED of F10214+4724 shows simul- 
taneously characteristics of almost edge-on and almost face-on emission. The failure of 
torus emission alone to fit the mid-infrared SED of F10214-I-4724 motivated Efstathiou 
(2006) to propose a three component model: an edge-on torus, clouds associated with 
the narrow-line region and a highly obscured starburst. The model assumed that the 
narrow line region clouds have a covering factor of 17% and two discrete values for the 
dust temperature (610 and 200K). The new MIPS and PACS data suggest that there is 
about a factor of two more power in the rest frame wavelength range 20-50/i?n compared 
with the prediction of the model presented in Efstathiou (2006). Considering that the 
new data samples the peak of the SED of F10214+4724 this implies that the covering 
factor needs to increase even more. It therefore seems likely that the magnification of the 
emission from the narrow-line region is considerably higher than that from the torus. We 
find that the new data can be fitted by assuming an additional component of narrow- line 
region dust at a temperature of 80K. The new fit is shown in Figure 1. The torus emission 
is assumed to be negligible in this fit. 
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